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Interesting and recommended websites: 

Lehrstuhl fur Didaktik der Mathematik Uni Erlangen, Prof. Thomas Weth 
www.didmath.ewf.uni-erlangen.de/Homepage/hp weth.htm 

Die Zeitschrift fur Mathematiklehrerlnnen (Universitat Salzburg und PH Salzburg) 
www.mathematikimunterricht.at/Newsletter/index news.html 


The Journal of Symbolic Geometry 
iournal.geometryexpressions.com/ 

Documents of Maths & Stats & OR 
www.ltsn.gla.ac.uk/headocs/ 

Among others a resource for background pictures 
staff.spd.dcu.ie./oreillym/geometry.htm 

Teach Engineering, Resources for K12 
www.teachengineering.com/index.php 

A publication of the Institute of Electrical and Electronics Engineers 
www.ieee.org/web/education/preuniversitv/tispt/lessons.html 

Online proceedings of the CADGME 2009 Conferrence is available at: 
www.risc.uni-linz.ac.at/about/conferences/cadgme2009/ 

The first Journal of Mathematical Modeling and Application of the Reference Center for 
Mathematical Modeling in Teaching - CREMM. The first number is on line at: 
furb.br/ois/index.php/modelling 

Teaching Math through Culture - a very interesting site! 
www.rpi.edu/~eglash/csdt.html 

Just recently I was informed by Philip Yorke (former DERIVE dealer for UK) about a 
very rich resource of A-Level Math (Exam questions and others). Many thanks Philip for 
this valuable notice. Have a look - or more than only one! 
www.mathsnetalevel.com 


Please inform me about interesting websites that we can share our resources, Josef. 
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Dear DUG Members, 


This is DNL#76 containing the remaining lectures of the special DERIVE session at 
ACA09 in Montreal. With DNL#77 we will proceed publishing original contributions 
intended for publication in the DERIVE Newsletter. 


The User Forum is not very extended in this issue because the articles needed a lot 
of pages. I received a mail from "Santa Claus" about a lot of DERIVE bugs. You will 
find them in the next DNL. 


As now WINDOWS 7 is on the market I have been asked several times if DERIVE 
will work under WINDOWS 7. Yes, it does. The respective information provided by 
DUG-Members is given in the User Forum. Many thanks to Gunter and Peter. 

Among many others Phil Todd joined the DUG this year. Phil is author of Geometry 
Expressions, which is an excellent piece of software, GE enables exporting of re- 
sults to DERIVE and to TT-Nspire as well. We will demonstrate this in 2010. 

We have good news for you. Bernhard Kutzler gave permission to put all Conference 
Proceedings which were published by bk-teachware on the ACDCA-website 
(www.acdca.ac.at). We will start uploading as soon as possible. Philip Yorke who pub- 
lished the Proceedings of the famous Krems Conferences (1992 and 1993) also per- 
mitted to upload the Proceedings. This will be more work, because we don't' have 
them in electronic form. I' II scan the books and produce pdf-files. (This will need 
some time, so please be patient.) 

Have a look at the recommended websites. I checked them just now - they all are 
valid at the moment. 


We wish a Merry Christmas and a Happy New Year 2010 

Noor and Josef Bohm 


Jose Luis Galan and Josef are planning 
TIME 2010 (right). 

Finally I'd like to remind you that dead- 
line for submitting a paper for the 
TIME2010 conference is 15 March. The 
website provides some information 
about accomodation and registration 
fees. 



The Conference website is www.time2010.uma.es. 


Download all ZWL-Derive- and Tl-files from 

http : / / www . austromath . at/ dug/ 
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Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. It 
must be said, though, that non-English 
articles will be warmly welcomed nonethe- 
less. Your contributions will be edited but 
not assessed. By submitting articles the 
author gives his consent for reprinting it in 
the DNL. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE & CAS-77 Newsletter will be. 


Next issue: March 2010 

Deadline 15 February 2010 

Preview: Contributions waiting to be published 

Some simulations of Random Experiments, J. Bohm, AUT, Lorenz Kopp, GER 

Wonderful World of Pedal Curves, J. Bohm 

Tools for 3D-Problems, P. Luke-Rosendahl, GER 

Financial Mathematics 4, M. R. Phillips 

Hill-Encription, J. Bohm 

Simulating a Graphing Calculator in DERIVE, J. Bohm 
Henon, Mira, Gumowski & Co, J. Bohm 

Do you know this? Cabri & CAS on PC and Handheld, W. Wegscheider, AUT 

Steiner Point, P. Luke-Rosendahl, GER 

Overcoming Branch & Bound by Simulation, J. Bohm, AUT 

Diophantine Polynomials, D. E. McDougall, Canada 

Graphics World, Currency Change, P. Charland, CAN 

Cubics, Quartics - interesting features, T. Koller & J. Bohm 

Logos of Companies as an Inspiration for Math Teaching 

Exciting Surfaces in the FAZ / Pierre Charland 's Graphics Gallery 

BooleanPlots.mth, P. Schofield, UK 

Old traditional examples for a CAS - what's new? J. Bohm, AUT 

Truth Tables on the Tl, M. R. Phillips 

Advanced Regression Routines for the TIs, M. R. Phillips 

Where oh Where is IT? (GPS with CAS), C. & P. Leinbach, USA 

Embroidery Patterns, H. Ludwig, GER 

Mandelbrot and Newton with DERIVE, Roman Hasek, CZ 

Snail-shells, Piotr Trebisz, GER 

A Conics-Explorer, J. Bohm, AUT 

Coding Theory for the Classroom?, J. Bohm, AUT 

Tutorials for the NSpireCAS, G. Herweyers, BEL 

Some Projects with Students, R. Schroder, GER 

Runge-Kutta Unvealed, J. Bohm, AUT 

The Horror Octahedron, W. Alvermann, GER 

RK6, Heinrich Ludwig, GER 

and others 

Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wiirmla, I ) I .list 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag. Josef Bohm 


The DERIVE-NEWSLETTER is the Bulle- 
tin of the DERIVE & CAS -77 User Group. 
It is published at least four times a year 
with a contents of 40 pages minimum. The 
goals of the DNL are to enable the ex- 
change of experiences made with DERIVE, 
TI- CAS and other CAS as well to create a 
group to discuss the possibilities of new 
methodical and didactical manners in 
teaching mathematics. 

Editor: Mag. Josef Bohm 
D'Lust 1, A-3042 Wiirmla 
Austria 

Phone: ++43-(0)6604070480 

e-mail: nojo.boehm@pgv.at 
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DERIVE and Windows 7 

Peter Hofbauer, Horn, Austria 

Hi Josef, 

DERIVE is running under Windows 7 without any problems. The Online Help cannot be shown (like 
under Vista). But there is a patch available: 

Download the respective patch from 

http : // www . microsoft . com/ downloads/ details . aspx?f amilyid=258AA5EC-E3D9- 
4228-8844-008E02B32A2C&displaylang=de 

(for 64-bit or 32-bit operation systems). 

Best regards 
Peter 

Note: Some screen resolutions may cause troubles (parts of the display are cut off and inserted at other 
places, . . .) Looks very strange. What to do? Decrease the resolution, then it should work properly. 

Gunter SchodL, Wr. Neusdtadt, Austria 

Hallo Josef! 

Derive runs properly under Windows 7. All what one has to do is installing a patch for the Online- 
Help. This can be found at www.microsoft.at under 

KB917607 

Windows6.l-KB917607-x64.msu for 64Bit Systems 
Windows6.l-KB917607-x86.msu for 32Bit Systems 

Best regards 
Gunter 


Gerhard Hagen, Friesach, Austria 


Dear all, 

I ask for help to solve a "DERIVE-mystery". I found the following function in an interesting article of 
the "Scientific News" and wanted to reproduce it with DERIVE: 



yj2a- x 

0<x<--T r 2 
2 R 


cl-T r 

— • 7^ _ x < v 0 • T r — — • T r 

Av(x) = < 

K 


pa ■ v 0 • T r - 2a ■ x 

v 0 -T R -~ T R 2 <x<v 0 -T R =d 


0 

d=v 0 -T R <x 


Av(x, a, tr, vO) := 

If 0 < x < a/2 * tr A 2 
7(2- a- x) 

If a/2 • tr A 2 < x < vO-tr - a/2-tr A 2 
#1: a-tr 

If vO*tr - a/2*tr A 2 < x < vO-tr 
V(2-a-v0-tr - 2-a-x) 

If vO-tr < x 
0 
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A TOOLBOX WITH DERIVE 

Alfonsa Garcia, Francisco Garcia 

Universidad Politecnica de Madrid. 

Departamento de Matematica Aplicada. E.U. Informatica 

alfonsa. garcia0eui. upm . es 
gmazario6eui . upm . es 

Gerardo Rodriguez 

Universidad de Salamanca. 

Departamento de Matematica Aplicada. E.P.S. de Zamora 
crerardo@usal . es 

Agustrn de la Villa 

Universidad Pontificia Comillas. 

Departamento de Matematica Aplicada y Computacion. ETSI (ICAI) 
Universidad Politecnica de Madrid 
Departamento de Matematica Aplicada. E.U.I.T. Industrial. 
avilla@upcomillas . es 


Abstract: An analysis of some of the characteristics of the European Higher Educa- 
tion Area (EHEA), its difference in relation to the current University teaching sys- 
tem and the role that new technologies might play in this new scenario has been 
performed. This paper suggests a new possibility in use of technologies: The de- 
sign of a “toolbox” with DERIVE instructions about topics in a usual Calculus 
course. 


Introduction 

The implementation of the EHEA (see [1, 2]) implies new teaching methods taking into account that 
the students now are the centre of the learning process. The role of teachers changes and they must be 
able to guide their students’ work (see [6, 7, 8, 9]). 

Teachers are currently being required to change the traditional teaching model in order to adapt to 
“learning based on competences”. It is necessary to define the competences to be acquired by the stu- 
dents after attending a course on a certain subject, and to design the activities according these compe- 
tences. 

Such change may fail if a considerable amount of effort, imagination, common sense, and hope is not 
devoted to it. The system inertia and the difficulty involved in designing effective activities must be 
borne in mind, since teachers generally have extensive experience in preparing expository lessons, 
with more or less of students, and although we have worked hard to find the best way of introducing 
and presenting concepts and results, we do not have any experience in guiding the search. 

The new teaching model implies more autonomous work by students. To see this it is merely neces- 
sary to analyze the new structure of studies, which are articulated in 60 yearly ECTS (European Credit 
Transfer System) credits, where each ECTS credit must reflect between 25 and 30 hours of the stu- 
dents’ “overall work”, only 10 to 15 of which must necessarily involve classroom attendance. 
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This is why teachers find it challenging to design a mathematical course for engineering students (for 
example a course on Calculus of One Variable or a course on Linear Algebra) taking these determin- 
ing factors into account. 

In our opinion we are doomed to design mathematical courses in which magisterial lessons (theory 
and problems solved in detail by the teacher), practical workshops (problem solved and Mathematical 
laboratories based in a CAS), and tutorial activities must be blended so that students can acquire the 
required competences. 

Laboratory classes must be designed with clear goals. Our proposal is that Computer Algebra System 
(CAS) could help in the automatic performance of certain tasks involved in the problem-solving proc- 
ess. For this, it would sometimes be necessary to use certain functions or commands which might al- 
ready be integrated to the system or which might have been prepared by the teacher, or even created 
by the students themselves. 

Working in a more autonomous way allows students to access computing technologies outside con- 
ventional training, which is why it is essential that they acquire the skills to make optimum use of 
them. The advantages of CAS must be boosted: visualization, computation facilities, the possibility of 
experimenting..., avoiding possibly damaging effects such as the lack of a critical attitude when con- 
sidering the computer response, an inability to interpret the results, etc. 

In any case, in order to be effective all the activities suggested must be designed without letting the 
intended goals out of sight, mainly taking into account the students at whom they are aimed. 

1. A box of mathematical tools 

Mathematical subjects, which are usually programmed within the first years of Engineering studies, 
have as their main goal the initiation of students into the language of Science and Technology and 
their preparation in the correct use of certain algorithms in problem solving. 

Teachers often complain that outside the context of the subject (for example, in later years) students 
are not able to use the mathematical skills acquired during the basic years of their training. They seem 
to have a kind of “mental laziness” that prevents them from remembering and using what they have 
learnt, and in many cases they do not have a fast and effective way to find the information or the ap- 
propriate methods either. 

The use of technology in the classroom is increasing and in certain cases students are even asked to 
define some tools that allow them to automate certain simple tasks, such as for example the implemen- 
tation of a function to calculate the tangent to a curve y — f (x) at a point and subsequently use it in 
other problems. In general, however, when teachers suggest to students that they should define a 
method aimed at automating a mathematical task they merely wish to help them understand the corre- 
sponding algorithm and they do not usually expect students to design their own resources or use the 
implemented functions in following years. 
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One educational activity is to encourage students to create their own well-organized “Toolbox” for 
solving mathematical problems. This toolbox is no more than a file or collection of files of utilities, 
programmed in the characteristic programming language of a CAS (Derive, Maple, Maxima) or even 
on a calculator with symbolic or graphic capacities such as the TI92 or Casio ClassPad 300. 

The teacher should suggest a series of essential tools, depending on the corresponding subject, that the 
box should contain. Students must to define the corresponding functions, test them, and add those they 
find appropriate. In addition, they must complete their work by writing a “brief user’s manual” for 
their tools. 

If they have created a good toolbox, and they have also understood the algorithms well, they will then 
have their own resource, which they will be able to use in other subjects during the same or following 
years. 

Before finishing with the general ideas on toolboxes it should be noted that nearly all CAS offer a very 
thorough toolbox and that the creation of new tools, or the modification of those already available, 
will only be necessary or convenient for reasons of ease of use or as a teaching strategy in order to 
adapt them to the user’s specific needs. 

2. A toolbox of Calculus 

In this section we propose a toolbox that can be done by students of Calculus in the first course of 
Engineering, whose reference text book may be [3 or 5]. 

Before present our toolbox, is necessary to make it clear that in DERIVE are implemented the most of 
the instructions for the study of Calculus of One Variable. Thus, the calculation of limits, derivatives, 
integrals or Taylor polynomials is “basic” using the DERIVE menu or instructions. 

Following the strategy referred to in the previous section, a supplementary toolbox has been created 
that can be used in a more specific way. 

With the toolbox our aim is to extend the use of the CAS, in our case DERIVE, and we therefore try to 
use, as far as possible, its symbolic, numerical, and graphic capacities. 

Some of the utilities are “improvements” to DERIVE instructions (for example the TANGENT2 method, 
where a distinction is made between differentiable and non- differentiable functions). Other tools are 
analogous to known instructions integrated in DERIVE, whose syntax is complicated. For this reason 
we prefer students to automate the algorithms according to their own criterion, so that they get a more 
continuous use when they have to apply these concepts in other topics. In addition, this automation 
“guarantees” that they have understood the corresponding concepts. 

Our toolbox could be designed with several compartments, one for each section of Calculus of One 
Variable course: 

• Complex Numbers. 

• Limits and Continuity. 
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• Differentiability . 

• Integral Calculus (including numerical integration). 

• Numerical Methods (for solving nonlinear equations). 

Below we briefly explain the tools. 


2.1. Complex Numbers 


The Complex Number compartment includes tools for: 

• Convert a complex number to exponential form. 

• Plot a complex number as a pair of real numbers. 

• Compute a list with the n th roots of a complex number. 

• Find the n vertices of a regular polygon. 

The following example allows one to calculate the n th roots of a complex number and, taking advan- 
tage of the graphic capacities of DERIVE, plot a start. 





x fin ^ A # 





> r 



raid -|*W=M*m\H’hW-l-|a|El«Ulil 




Figure 1 : A regular hexagon with a six-pointed start 


2.2. Differential calculus 

The utilities performed allow the analysis of the continuity or discontinuity of the function at a point 
(studying the one sided limits, analyzing the equality of the obtained values and comparing with the 
value of the function at the point). The differentiability of a function at a point is studied following the 
same strategy. 

The calculation of the tangent line has been implemented, with additional information to the TANGENT 
function of DERIVE, since it tells us when the function is not differentiable at the point (see figure 2). 
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The tangent 

#1: f(x) := 

f(a + h) - f (a) 

#2: DL(a) := lim 

h->0- h 

f(a + h) - f(a) 

#3: DR (a) := lim 

h-*0+ h 

TANGENTE2(a) := 

If DR (a) - DL(a) = 0 a (l/DR(a)) A 2 > 0 

#4: y = f(a) + DR(a)-(x - a) 

"f has no derivative' 1 

f(>0 := 

If x < 0 

#5: x 

x a 2 

#6: ilMdaHiiMiM 

#7: f has no derivative 

With the DERIVE instruction TANGENT the result is not correct 
#8: TANGENT (f(x), x, 0) 

#9: x 


Figure 2: TANGENT and TANGENT2 functions 

Taking advantage of DERIVE graphical capacities and of the structure of the IF instruction, the 
monotonicity and concavity or convexity of a sufficiently differentiable function can be analyzed. 

For example, with the instruction CRECE(x):= IF(F’ (x)>0,F(x)) it is possible to represent the 
curve y = F(x) in the intervals where F is increasing. 

The graph of figure 3 has been obtained by applying the instruction CRECE(x) with F(x) = x 3 - x 



Figure 3: Intervals where F(x) = X 3 — X is increasing 
Other functions about Differential Calculus in the toolbox may be: 
R0LLE_P0INT for finding a point according Rolle’s Theorem hypotheses. 


LAGRANGE_POINT for finding a point according Lagrange’s Theorem hypotheses. 


II 
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These functions can work in exact or approximate way, because it is necessary solving equations. 

2.3. Integral calculus 

In order to introduce the Riemann integral we can use (see [4]) the DERIVE instructions: Areallnder- 
Curve and LEFT_RIEMANN. Furthermore we have included in the toolbox the DER_RIEMANN for calcu- 
lating the sum for right rectangles associated to Riemann sums. We also have implemented the 
rect_izq and rect_der functions for plotting the left and right rectangles associated to Riemann 
sums 

We promote our students to define simple instructions for computing lengths, areas and volumes with 
a syntax more recognizable that POLAR_ARC_LENGTH , PARA_ARC_LENGTH , POLAR_AREA, 
V0LUME_0F_REV0LUTI0N , AREA_0F_REV0LUTI0N , etc. 

For instance, to calculate the length of an arc of the curve y — f (x), it is possible to define the func- 
tion LEXP and theoretically find the length of any arc of curve. 

The figure 4 shows the DERIVE implementation and the calculation of the parabola’s length y = x 2 
between the abscissas 1 and 2. 


#1: 

fO) &- 




b 







#2: 

LEXP (a, b) := J + f'O) ) dx 




a 



The 

parabola length 




2 



#3: 

fO) := x 



#4: 

LEXP (1, 2) 





LNC/S5 - 2-07 + 4-0 - S) 

ft 

#5: 


+ 0.7 

— 



4 

2 

And 

the approximate value 



#6: 





Figure 4: The length of an arc of curve 

Finally, students can define tools for numerical integration using the Composite-Trapezoidal rule and 
Composite-Simpson rule (see figure 5). 
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# 1 : 


#2: TRAP (a, b, n) := 


b - a ( f(a) + f(b) n - 1 


A( 




#3: SIMPSON(a, b, n) := 


-a f n/2 f 

If (a) + f(b) + 4. Z fl a 


(2-i - 1) • (b - a) 




#4: 


fGO 


- 2.X 

x -e 


2 

x +4 


#5: 


1 


J 

-1 


f (x ) dx 


# 6 : 


1 


- 2-x 

x -e 


2 

x +4 


dx 


#7: [TRAPC-l, 1, 10), SIMPS0N(-1, 1, 10)] 

# 8 : 


[-0.4342879366, -0.4215187544] 


Figure 5: Numerical Integration 


2.4. Numerical calculus for solving nonlinear equations 


The methods usually taught to students for solving equations are: Bisection, Newton and fixed-point. 
DERIVE has the NEWTON and FIXED_POINT algorithms integrated, so we can propose that students 
define a function to implement the bisection method. An algorithm for this method could as follows: 


H(a, b) := 

If F((a + b)/2)-F(a) < 0 
[a, (a + b)/2] 

If F((a + b)/2)-F(b) < 0 
[(a + b)/2 , b] 

[Ca + b)/2, (a + b)/2] 

BISECC(a, b, n) := ITERATES (H(v , v ), v, [a, b], n) 

1 2 


Figure 6: Bisection Method 

The instruction BISSEC(a,b,n) (see figure 6) computes the interval obtained after n iterations of Bi- 
section method, to solve the equation f(x) = 0 in the interval \a,b ] . 


3. Using the toolbox to solve technical problems 


Students are encouraged to use the toolbox to solve technical problems. Here we set out one example, 
taken from an Environmental Sciences exam that has been used by our engineering students. 

Specifically, the equation to be solved (which appears together with its solution in the figure 6) pro- 
vides the height of the chimney of a thermoelectrical station in a rural area with the required quality 
standards. 
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Numerical methods 

#1: FOO := 

#2: PUNTO_FIJO(xO , n) := ITERATES(F(>0 , x, xO, n) 

#3: FOO := 


#4: 


NWfa, n) := ITERATES 


F&O 

. x . a. n 

F'&O 


H(a, b) := 

If F((a + b)/2)-F(a) < 0 
[a, (a + b)/2] 

#5: If F( (a + b)/2)-F(b) < 0 

[(a + b)/2, b] 

[(a + b)/2 , (a + b)/2] 


BISECCfa, b, n) := ITERATE5(H(v , v ), v, [a, b], n) 
# 6 : 12 


Taken from an exan in iiy University (April 200Q) . 
Topic: Environmental techniques 


The goal: To find the iiiniiiuii height of a chimney of a thermoelectric central 
in a rural area to be concordant with the quality environmental standards. 


The equation to be solved for finding the height: 
H=l. 6*575. 29 a (1/3)(1 0(S88- H) a (2/3)|)/5 

1/3 2/3 

1.6*575.29 *(10.(588 - x)) 

#7: F(x) := x 

5 


1/3 2/3 

1.6-575.29 .(10.(588 - H)) 

#8: F(H) := H 

5 



Sketch the graph to approach the initial values of bisection method: 
#9: BISECC(200, 500, 10) 


# 10 : 


200 

500 

350 

500 

350 

425 

387.5 

425 

387.5 

406.25 

396.875 

406.25 

401.5625 

406.25 

401.5625 

403.90625 

401.5625 

402.734375 

402.1484375 

402.734375 

402.1484375 

402.4414062 


Figure 7 : The height of a chimney 


For engineering purposes the chimney height is around 400 meters. 
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Conclusion 

The use of a toolbox adapted to each student’s individual needs containing instructions that the stu- 
dents themselves find useful must, we believe, surely reinforce the learning process, since students 
participate actively and the use of the CAS is not limited to what is sometimes “blind” use of the CAS 
instructions. 
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An Implementation of the Mayan 
Numbering System in DERIVE 

Eugenio Roanes Lozano 
Francisco A. Gonzalez Redondo 

Algebra Dept., School of Education 
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Abstract 


The Mayan number system is a base 20, positional (to be read from top to bottom, not 
from left to right) system that makes use of a symbol representing zero. It has slightly dif- 
ferent variations when used for counting days (in religious and astronomical contexts). 
Therefore, 20 symbols are needed to represent 0,1, 2,..., 19. Of these, the zero was denoted 
by a shell and the positive ones were represented using dots and horizontal segments. If a 
number is greater than 20, the symbols corresponding to units, twentieths, 400’ s, 
8000’ s... are stacked from bottom (units) to top in pure base 20. We have implemented a 
procedure that allows to convert numbers between any bases, and that returns the output 
in (row) vector style; another procedure that builds the 20 Mayan symbols for 0,1, 2,... 19; 
and yet another procedure, that uses the previously mentioned procedures, converts any 
number from base 10 to base 20 and represents it in the Mayan numbering system. We 
believe this is an interesting example of synergy among different branches of knowledge 
(Mathematics, History of Mathematics and Computer Science), that can increase the in- 
terest of the students for different topics. 


1. Introduction 

Let us try to briefly describe the evolution of the Spanish educational system at Secondary 
Education as regarding Mathematics. Until the mid 70’ s, the curriculum focused on classic 
plane geometry. After the arrival of Modern Mathematics, set theory, relations, correspon- 
dences, algebraic structures, number base changes, etc. were studied instead. Now the curricu- 
lum of 10-16 year old students focuses on practice and experimentation, oriented to achieve 
basic skills. 

Nowadays, possibly as a consequence of an incorrect implementation of the latter ideas, 
many freshmen at university do not know or do not master basic concepts, procedures and 
algorithms. Some Schools, like the School of Mathematics of our University, have introduced 
“0 courses” devoted to those students that can’t follow the subjects of the first year at univer- 
sity. 

For example, they can’t perform number base changes, because they have never studied 
this topic. 
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We propose here a unit at Teacher Training level organized in three sessions of 1.5 hours 
about number base changes. Here, the Mayan numbering system is introduced as a historic 
justification for the need of working in other bases. Another justification is the use of base 2 
(and, indirectly, base 16) by computers. Therefore, this is a multidisciplinary unit, involving 
Maths, History and Computer Science. The working scheme for the unit is the following: 


Studying the problem (including the historic antecedents) 

i 

Translating the problem into mathematical language (mathematization) 

I 

Solving the problem 

i 

Implementing the solution found 


There is a didactic antecedent to this unit, at a lower educational level and without using 
technology: (Roanes, 1971). We do not know of any other similar implementations in any 
Computer Algebra System (CAS), apart from (URL, 2009), that uses special facilities for in- 
serting captured drawings and doesn’t have a didactic orientation, or the very similar (Roanes 
Lozano & Gonzalez Redondo, 2009), by the same authors of this paper, but written in the 
CAS Maple (that was successfully experimented during the 2008-2009 course). 


2. The Mayan Numbering System 

The Mayan civilization developed a characteristic world of astronomical and mathemati- 
cal knowledge in (and from) their temples, the most appropriated place for a political system 
organized as priestly autocracy which based its power upon Astrology. 

Without any observational instrument (they did not know glass and, consequently, could 
not manufacture lenses), without sand or water clocks for computing periods of time as hours, 
minutes or seconds, their smallest unit of time was the kin or day, whose duration was meas- 
ured through the repeated registration of the shadows projected by a wooden bar placed verti- 
cally on the ground, known as the gnomon (also used by the Greeks for determining the azi- 
muth and the height of the sun). 

Taking the day as basis, the subsequent units of time where: 

• the uinal or month of 20 kins (days), 

• the tun or year of 1 8 uinals (18-20 = 360 days), 

• the katun or period of 20 tuns ( 20 • 360 = 7200 days), 

• the baktun or period of 20 katuns ( 20 • 20 • 360 = 144.000 days), 

• the pictun or period of 20 baktuns ( 20 • 20 • 20 • 360 = 2.880.000 days), 

• etc. 

As it seems obvious, Mayan astronomers knew the corrections that had to be made after 
considering a year of 360 days instead of the solar year of 365.25 days (Grube, 2000). 
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With this system, and counting long periods of time, they estimated the phases of Venus 
(which they rounded up to 584 days). Mercury and Mars; the duration of a lunar month in 
what today we would write as 29.53020 days; and they even registered periods up to 
300.000.000 years. We can confirm all these through their commemorative steles and the 
pages of the few Mayan books ( codices ) which have being preserved (Landa, 1937), where 
we find an original numbering system with three basic features: it had a base 20, used the 
place-value system and introduced the zero. 

In fact, in association to its religious origin, numerical quantities registered by the Mayas 
where strictly joined to the god corresponding to each order: the god carrying kins, the god 
carrying uinals, the god carrying tuns, etc. Their conception of the “zero” was based upon this 
religious nature: it was a symbol created in order to fill any possible gap in the place where 
some numerical quantity should be supported by the god corresponding to the said order 
(Ifrah, 1999). 

In short, this system conceived at the temples was perfect for representing dates (elapsed 
days), prepare ritual celebrations, etc., although with it no arithmetical operation could be 
performed. For such operations which were well apart from religious rites and more related to 
administration and commerce, the Mayas put aside the irregular tun (year) of 360 kins, assum- 
ing the natural order 20 • 20 = 400 proper of a base 20 numbering system. 

While in our usual base 10 system we need nine numerals (and the zero), for writing their 
numerals the Mayas needed nineteen numerals. But the symbols used for each of these nine- 
teen units were very simple: dots and horizontal lines (see Figure 2). 



Figure 2: Mayan numerals. 


When we write today a positive four-numerals integer n in our usual base 10 system, n = 
wxyz, where vv, x, y, y z can vary from 1 to 9 (or 0), we are really shortening the expression: 

n= w-10 3 +X-10 2 + y -10 1 + z -10° 

For writing any number greater than nineteen, the Mayas placed the symbols referred in a 
vertical column containing as many rows as orders of units, including a very special symbol 
(reminding us of a shell) for filling up an order that could not be left empty of dots or lines: 


[i] 


The pictures were taken at my 3 weeks travel through Mexico in November 2009, Josef. 
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the Mayan zero. Through this notation the Mayas could write as follows what in our sys- 
tem would be a number m : 

m = ... + w • 20 3 + x ■ 20 2 + y ■ 20 1 + z ■ 20° 

For example, the number: 

1 , 368, 280 = 8 - 20 4 + 11 - 20 3 + 0 - 20 2 + 14 - 20 1 +0 - 20 ° 

was written using dots, lines and shells as it is shown in Figure 3. 



Figure 3: Example of Mayan numerals. 

For those interested in further mathematical developments of the Mayas see the book by 
Romero (2004). 


3. What DERIVE does about number base changes 

Regarding number base changes, the CAS DERIVE only provides the possibility to choose 
the input and output bases (from 2,8, 10 and 16) in: 

Options > Mode Settings 

Therefore we’ll try to complete its possibilities implementing the corresponding proce- 
dures so that it can convert any given integer number into Mayan notation. The mathematical 
background of each procedure is previously taught or retraced, and the students are guided 
through the development of the implementations (these details are omitted here for the sake of 
brevity). As said above, this unit was successfully experimented during the 2008-2009 course 
with the CAS Maple. 
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4. A DERIVE implementation for number base changes 

It is straightforward to implement the auxiliary procedure (or function) “integer quotient” 
in DERIVE. As DERIVE’ s command MOD returns the remainder of the integer division, we 
can follow the integer division definition: 

a - M0D(a, b) 

quo(a, b) := 

b 

Now the previous procedure can be used as subprocedure by a short recursive procedure 
that converts numbers from base 10 to any base, and returns the output in (row) vector style: 

base_bl(n, b) 

If n < b 
[n] 

APPEND(base_b1 (quo(n, b) , b) , [M0D(n, b)]) 

Examples: 


base_bl (399 , 20) = [19, 19] 
base.bl (400, 20) = [1, 0, 0] 


5. A DERIVE implementation of the Mayan numerals 

Representing the Mayan numerals (0 to 19) is tricky and the less interesting procedure 
from the mathematical point of view. 

mayal9(n) := 

Prog 

If n = 0 

DISPLAY(" 0 ") 

If n = 1 

DISPLAY(" . ") 

If n = 2 

DISPLAYC" . . ") 

If n = 3 

DISPLAYC" ... ") 

If n = 4 

DISPLAYC" ") 

If n = 5 

DISPLAYC" ") 

mayal9_auxCn) 


RETURN 


ll 
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mayal9_aux(n) 

Prog 

If MOD(n , 5) > 0 
mayal9(M0D(n , 5)) 

If quo(n, 5) = 1 
mayal9(5) 

If quo(n, 5) = 2 
maya_10 

If quo(n, 5) = 3 
maya_15 
RETURN " " 

maya_10 :z PROG(DISPLAY( ), DISPLAYC )) 

maya_15 :z PROG(DISPLAY( ), DISPLAYC ), 

DISPLAYC )) 


Example: 


* 9 : 


mayal9(ll) 


6. A DERIVE implementation of a Mayan number system converter 


Once we have the procedures of Sections 4 and 5, it is very easy to implement a Mayan 
number system converter. The idea is to convert first the given number to base 20, to trans- 
form the numerals to Mayan notation, and, finally, to stack these numerals. 

mayaCn) := 

Prog 

maya_auxCbase_b1 Cn , 20)) 

RETURN " " 

maya_auxCR) 

Prog 

If R t [] 
mayal9CRrl) 

If R * [] 

DISPLAYC" ") 

If R t [] 

maya_auxCRESTCR)) 

RETURN " " 
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Example: 2721 = 6 x 20 x 20 + 16 x 20 + 1 


#10: maya(2721) 
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8. Conclusions 

We do not know of any other similar implementations in any CAS, apart from (URL, 
2009), that uses special facilities for inserting graphics, and the very similar (Roanes-Lozano 
& Gonzalez-Redondo, 2009), written in Maple. 

We believe this is an interesting example of synergy among different branches of knowl- 
edge (Mathematics, History of Mathematics and Computer Science), that can increase the 
interest of students for different topics. 

The only negative issue we find is that parts of the implementation are a bit tricky (Section 
5), something that can distract students from the mathematical goals of the unit. Perhaps this 
part could be provided to the students implemented beforehand. 



The "Observatory", Chichen Itza 


Uxmal, Temple of the Wizard 
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Abstract 

Even in times of a CAS a lot of basic manipulating skills - up to a not commonly 
agreed amount - seem to remain necessary in mathematics education. The 
teacher can provide programs and functions to encourage self responsible learn- 
ing, training and repeating these skills. It is the teacher - and sometimes the cur- 
riculum and the given "standards" - which level of difficulty is appropriate. 

I will show some examples reaching from calculating GCD and LCM via factoriz- 
ing and expanding expressions up to basics of calculus using various computer 
algebra systems. This could raise a discussion about possible fields in math cur- 
riculum which might be suitable for CAS-supported exercising. And this could 
also raise a discussion about advantages and disadvantages of using these or 
similar tools. 


Let me start with a view "some years" back. It was in 1994 when I received a diskette from 
my friend Jan Vermeylen (from Belgium) - no email at these times- containing three files 1101 : 

QUADR_EQ . MTH 

POLY_FAC . MTH 

and 


REKRIJ.MTH 

■ ■ > 

jj] QUADR_EQ 1 KB r4TH-Datei 

■j REKRIJ 2 KB NTTH-Datei 


16 . 07 . 199415:02 
14 . 07.2003 10:35 


Let's try if REKRIJ works with the recent version of DERIVE, too: 
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COMMAND: T'H.W i'TTSl Build Calculus Declare Expand Factor Help Jump soLue Manage 

Options Plot Quit Remove Simplify Transfer Unremove moUe Window approX 
Enter option 

Sinp<ttl9> D:\DOKUS\SCHULE\SKI Free:100x Derive Algebra 


#19: SER(5) 


" 15 

de element is 

. 15 

de som is : 

" 12 

de element is 

. 5 

de som is : 

" 4 

de element is 

. 19 

de element is 

" 6 

de element is 

. 8 

de som is : 

" 7 

de element is 

. 9 

de som is : 


270 1 " 

2055 J 

30 1 

60 J 
26 1 
146 J 

12 1 
48 J 

117 1 
747 J . 


REKRIJ with DERIVE 3.14 and with DERIVE 6.1 


I left the file in its original Flemish (= Dutch) version. All problems are connected with arith- 
metic series. Examples like these could and still can be found in nearly all textbooks. 

The first example should be read as: the 15 th element of an arithmetic series is 270 and the 
sum of the first 15 elements is 2055. Find the series (1 st element and difference). 

Solutions were not provided. I will come back to a modern version of this tool later in the pa- 
per. 

There was another attempt providing training tools for elementary calculation skills presented 
by Heinz Rainer Geyer [12] : 


SQUARES (5, 20 > 



1138 : MIXED CIO, lOOJ 


1139 : CHECK 


1634 


38 
43 = 48 

34 = 34 
103 = 103 

4 = 4 

35 = 35 
3640 = 3640 

-61 = -61 
-4 = -4 
139 = 139 
L -46 = -46 


48 , 


1634 

38 

1224 

36 

37 ,p + " 

66 

59 ,p - ,p 

55 

16 ” + ” 

19 

52 

70 

36 ,p - ,p 

97 

66 ,p -' p 

70 

43 ” + ” 

96 

42 

88 

1224 

34, 


=34, 37 +66 = 103, 59 - 55 = 4 , 16 + 19 = 35 , 


"wrong'' 
"r ifirlrt'" 
,p r isfli-t" 
,F r ight" 
"r ight" 
,F r ight" 
,F r ight ' 1 
”r isfli-t” 
”r ight ' 1 
"r ight" 


Sometimes I like to reanimate my old DOS-versions of DERIVE. 
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Allow a time step to the handheld devices: 

My school was the first College for Business Administration in Austria where CAS was used 
in all classes. When the TI-92 appeared on the market we - the math teachers - asked the 
parents to buy a device for their children for enabling a modern and attractive math educa- 
tion. It was in these times when I remembered another old BASIC-program which was later 
replaced by a DERIVE tool for exercising solving quadratic equations applying the rule of 
Vieta' 231 . 


2 

vi'etaO = ( z - 14 -z = -13) 
ant(l, 13) = ri'chti'g! 

2 

vi'eta() = ( z + B-z = 0) 
ant(0, -B) = ri'chti'g! 

2 

vi'etaO = O' - 15-n = -48) 

ant(-3, -16) = falsch, sondern: 3, 16 


2 


vi'etaO = (y 
ant(S, -9) = 

results 0 = 


+ v = 72) 
ri'chti'g ! 

Aufgaben gestellt: 
davon ri'chti'g: 
das sind in %: 


4 " 
3 

7B _ 


And here is the Tl-Version: 


|-: J ; ,|| l |'| 

■■ ' 

' ! ! v '' 


' #r ' 




s I c 


P:- 1 0 

: .S L =!■■! :i ' I! .. 


Problems: 1, correct 1 

g^ + 18 ■ g + 80 = Q 

1. Solution: 

-10 

2. Solution: 

-8 

right 

End = ESC. next = any 

MAIM DEG AUTO FUMC 0/30 ~ um \ 


fVj 1 ; 1 1 1 r (i 

■■ ' 

' ! ! v ' 


' #r ' 




0 5 ! c 


P:- 1 0 

: .S L =!■■! :i ' I! .. 


Problems: 3, correct 3 

C 2 - 10 c - 11 = 0 

1. Solution: 

-11 

2. Solution: 

1 

sorry, false xl = 11, x2 = - 1 
End = ESC. next = any 

MAIM DEG AUTO FUMC 0/30 ~ IJIHMl 


The idea remained the same, the tools changed. The pupils used the program during the 
breaks in school, often in the bus or railway on their way to or from school and sometimes 
only just for fun. One female pupil told me that she used her "VIETA" to concentrate before 
learning for other subjects. 

Just to demonstrate that the basic ideas can be realized using (almost) all tool I tried to pro- 
gram my “VIETA” with WIRIS, a CAS which recently has become popular in some Austrian 
schools. It is not the tool which decides what can be done, it is always the user (= teacher) 
who forms the software according to his or her didactical intentions and needs. 

See the WIRIS screen shot on the next page. 

I'd like to know in how many CA-systems I will program VIETA in the future? 
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KIM WIRIS Desktop 2 - D:\DOKUS\wiris\vieta.wiris 


B WIRIS CAS 


File Edit Tools Help 


Edit Operations Symbols Analysis 

Matrix Units Combinatorics Geometry Greek Prograr 


El 

% X taj 

Ak 


ip 

3 % 


1 




n 


Let's have another jump: Even in times of 
CAS the basic rules of differentiation and 
integration should be mastered by the stu- 
dents. Calc() - a program package for the 
Tl-handheld devices CAN help 111 . 

It is in your responsibility as teacher up to 
which level you think that the manipulating 
skills in calculus are indispensable. 


F2- 

F> 

FH- 

F£t 

Dihh 1 

Dihh II 

IntiSK-il 1 

InUdKOl II 


calcO 




DIFFINT 

ftflD EI-IhlT 

FUNC 

1/30 




i if i[inUJKoi u[ 

■*gr n. t f?- a 

f t A 1 ft 1 fi.ip' 1 1 1 1 ft frifi- 1 fr* fill 1 1 Ti I 1 £- y li "1 111 


Virr ■lIlrM.'s Jr y ■ HI I r| Jr -y 1 ■ HHl 1 


2:1 Power functions 
3s Product Rule 
4 : Quot i ent Ru 1 e 
5: Chain Rule 1 
6: Chain Rule 2 
7- Prod* Ri & Ch, R* 
SlQuot, Ri & Ch* R* 

2! i Power f unc^ons 
3! Product Rule 
4: Quotient Rule 
5: Chain Rule 
6:Diff Quiz 



_ _ unctions 

3 s Subst i tut i on 
4: Int. bu parts 
5 s Part i a I F ract ions 
6 : Integration Quiz 


F5- 

ihU*r«i u 


l: Polynonials 
2 : 2 Po war f unct i ons 
3 * Substitution 
4s X by parts/uo.S. 
5:X by partsAj-S, 

6 : Part i a 1 Fract ions 

iMfflBMMa 


All tasks are random generated. I'll demonstrate the basic idea by using the respective 
DERIVE tool which was created when we had the opportunity to teach math in the PC lab 
working with DERIVE. 
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I provided one program for training the differentiation rules and another one for exercising the 
integration rules. Both of them are offering two levels of difficulty. This is the "start"-screen for 
integration. 


2 lev 

els 

of diff 

iculty 

: (1) 

and 

(2): 



Level 

1: 








poll= 


, polynomials. 



pfl= ... 

I of power 

functions 

subst 

1= 

. . . Int. 

by su 

bst-itL 

jtion 

, partl= . 

. . Int. by 

sarts 

pfrl= 


, partia 

1 frac 

tions 


quizl= . 

, . random p 

roblem 

n tas 

ks 

of same 

type: 

pollb( 

In). 

pflb(n), .. 

, , quizlb(n 


Level 

2: 








pol2= 


, polynomials, 



P f2= ... 

I of power 

functions 

subst 

2= 

. . . Int. 

by su 

bst-itL 

Jtion 

, part2= . 

. . Int. by 

parts 

parts 

= ■ 

. . Int. 

parts 

& inc, 

sub 

stitution 



pf r2= 


, partia 

1 frac 

tnons , 


quiz2= . 

, . random p 

roblem 

quiz= 


, random 

probl 

em out 

: of 

all levels 



n tas 

ks 

of same 

type: 

pol2b( 

^n) , 

pf2b(n) , .. 

, , quiz2b(n 

) 

check 

(Ir 

itegral)= 

or ch 

eck([] 

inti ; 

Int2 ; ...] 

) gives the 

answer (correct/true), 

ans= 

ret 

:urns the 

cor re 

ct antider 

ivative 




Let's have 5 integrals of level 1 to exercise integration by substitution together with the cor- 
rect answers. The students should do the integration by hands and then double check their 
results. 



The next figure will show a level 2 quiz (5 problems). This is a random selection of all possi- 
ble problems. 

When I was a teacher and needed problems for a test I did not prepare the problems at 
home but I ran my problem generator files and eg simplified twice quizlb(IO) for two groups 
of students in the class. So they could not blame me for providing extra difficult tasks for the 
test - the computer made the selection. 
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1 

lS-Jx + 

2 

4-x 


12.x 


5 


#54: quiz2b(5) = 



54 


5 3 2 

56 -x 4 x 7-x 5-x 1 

5 5 5 4 3 


#56: check 


6 

56- x 


6 


3/2 

15.x -2 1 



x 

3 


= wrong 


Do you find the mistake? 


DERIVE has a wonderful feature - its "STEPWISE SIMPLIFICATION”. Unfortunately no 
other CAS which is on the market now has a similar option. If the student has a wrong inte- 
gral and he/she does not know how to obtain the right answer then he/she can ask DERIVE 
to stepwise integrate. 


I'll illustrate this applying stepwise integration on 
the second problem of the list of integrals by 
substitution from above: 


#50: 


4-x - 2-x 

20 -e ■ (1 - x) dx 


J a-F(x) dx -> a-J F(x) dx 


#51: 


- 2-x + 4-x 


20 - 


■ (x - 1) dx 


a-F(x) 

a ■ F (x ) e 

e ■ F r (x) dx — > 


2 

#52: - 2-x + 4-x 

5-e 


#53 : 



There is a similar DERIVE program for exercising, repeating or deepening the differentiation 
rules. Most of my programs which were produced in the last years show the same structure. 
The user hast to load the program as a utility file and then simplify the command start. Then 
he is presented the instructions how to run the tool. 

It might be a good moment now to stop demonstrating more examples but asking which 
fields in school mathematics might need some exercising? 

Which manipulating (and other?) skills (techniques, competences?) could (could not), should 
(should not) be trained, deepened and/or extended supported by CAS (or other technolo- 
gies)? 

My list is more or less CAS-oriented but I hope that I can extend it to other technologies and 
that I will be able to reach one or the other goals of higher competences than "only" manipu- 
lating skills. 
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This is my list: 


Working with fractions 

Expanding and Factorizing Expressions 

Power Rules 
Quadratic Equation 

Completing the Square 

Simultaneous Linear Equations 
Working with complex numbers 

Long Division of Polynomials 
GCD & LCM 

Working with units (time, length, area, ....) 

Basic Problems with Linear Functions 

Working with Vectors and Matrices 

Set Theory 

Truth Tables (Boolean Expressions) 

Analysing Conics 

Basic Problems for Financial Mathematics 

Solving Triangles 

Limits 

Arithmetic & Geometric Sequencies and Series 

Do you remember? Jan Vermeylen 1990! 

Investigation of Sequences 
Differentiation and Integration Rules 

Discussion of Curves 
Taylor Series, Fourier Series 
Typical Forms of Differential Equations 
Implicit Differentiation 
Conversion between number bases 


Recognising Function Types from their graphs 

Sketching derivative and/or antiderivative to a given function graph 


Finding Polynomial Functions 


All the items in bold letters are - more or less - ready. The underlined ones do not provide 
randomly generated problems but help solving given problems. All the items in italics are on 
my TO-DO-list. 

I in the following I will show a selection of screen shots. All tools from above are done using 
my favourite CAS, DERIVE. But it is possible to convert all ideas to any other CAS which is 
programmable like MATHEMATICA, MAPLE, MAXIMA, WIRIS, TI-NspireCAS, ... 

In my abstract I mentioned GCD and LCM. So I'll start with exercising calculating GCD and 
LCM of numbers: 
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#11: start 

Trainer fur ggT kgV und kgV - Trainer for gcd and 1cm 


ggtkgvz liefert eine Aufgabe - gcdlcmz offers one task 


res= gibt die richtige Antwort aus - res= delivers the correct answer. 


Find gcd and 1cm of 

#12: gcdlcm = 

. [36323, 463, 1300, 13972] _ 

GCD = 36 

#13 : res = 

_ LCM = 406949400 . 

See next the the "modern" version of the arithmetic series program (remember Jan Vermey- 
len's file from the DOS times!) This version provides the solutions, too. 



In 2000 bk-teachware published my Mathe-Trainer containing packages for the TI-92. There 
was one package algebra() which was intended for my students to exercise - among others - 
expanding powers of binomials 111 . Some years later when I could use DERIVE, version 5 was 
on the market and with it the possibility for programming. 

Two utility files were produced: 

expand. mth and factor.mth 
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#1: L0AD(D: \DOKUS\SCHULE\ski ns\expand .mth^J 

#2: start 

Expanding Expressions 


sq= gives a task for squaring a binomial 

cu= gives a task for cubing a binomial 

tr= gives a task for squaring a trinomial 

pr= asks for a product of two binomials 

di= sum times difference of two monomials 

quiz= presents a random problem of the above tasks 


sqb(n), cub(n) quizb(n) gives n tasks of the same type 


res?= gives the correct answer and ch(answer)= analyses single results 
The students are offered the instructions. 

This might by the start of a session exercising cubing binomials: 

3 

#4: cu = - (2-n + 3-u) 


3 2 2 

#B : ch(- S-n - 36- n ■ u - B4- n- u - 9) = check cubes! 

3 2 2 3 

#6: res? = - B-n -36-n-u-B4-n-u -27-u 


#7: cub(3) 


3 

(P + 7-s) 


3 

(2-g + v) 


3 

(S-m - t) 


3 2 2 3 

p + 21- p -s + 147- p-s + 343- s 

3 2 2 3 


#8: res? = 


S-g +12-g -v+6-g-v +v 


3 2 2 3 

. B12 - m - 192 -m -t + 24-m-t - t 

As you can see, I tried to include some error analysis for the most common types of errors. 

(This is a nice programming challenge - specially collecting the "common types of errors" 
together with the students.) 
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#1: LOAD(D:\DOKUS\SCHULE\skn'lls\factBr.mth) 

#2: start 

Program for training factorizing expressions 


II 

Q 

4- 

, . factor out, 

P f= .. 

. partially factor out 

sq= . , 

, . complete square 

cu= . . 

. complete cubic 

di= . , 

, . difference of squares 

cc= , . 

. sum or difference of cubics 

tr= . . 

, . trinomial 

qu iz= 

. . . random problem 


With fob(n) , pfb(n) , sqb(n) quizb(n) n tasks of the same type will be offered 

ch= shows the correct result 


poly= gives a factorizable polynomial 

polyb(rO= gives n polynomials 

ch= factorizes for rational zeros, 

chr= for irrational and chc= for complex zeros 


#3: 


quizb(B)= 


2 10 2 2 7 5 2 4 8 

-630- p -s -t -1680- p -s -t -U_20-p -s -t 

624 B 2 B 426 

-360-h -u -w-180-h -u -w+B40-h -u -w 

3 2 4 2 B 

432-f -k-1 +768- f ■ k -1 +336-f-k-l 

333 263 93 

67B-c -1 -n +13B0-C -1 -n -t67B-c-l -n 

4 2 3 2 4 3 

80- b ■ u -y +B00- b ■ u -y 

2 4 2 3 3 2 1 


-70-p -s -t -(3-s -+4-t ) 


#4: 


ch= 


4 2 4 

180-h -u -w ■ (w-h) ■ (2- h+3- w) 

2 2 3 3 

48-f-k-l ■ (9- f +16- k +7-1 ) 

3 3 3 2 

67B-C-1 -n -Cc+1 ) 

2 2 3 2 2 

20-b -u -y * (4- b +2B-u ) 


The instruction informs that this file can also be used for training factorizing polynomials on 
different levels of knowledge: starting with only rational zeros one can proceed to irrational 
zeros and then finish with complex zeros depending on the age of the students. 
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#9: polyb(B): 


B 4 3 

x -B-x +S-x 


3 2 

x ■ (x -S ■ x+8) 

6 B 4 3 

-x -lB-x -4B-x +64 -x 


3 2 

x ■ (1— x) ■ (x+B) 

6 B 4 3 

36 -x -2BB-x +234 -x +10B-X 

#10: ch= 

3 

6-x ■ (x-6)- (2-x-3)- (3-x+l) 

3 2 

B-x 104 -x 194 -x 2 OB 


2 

2- (B-B-x)- (4-x -4-X+13) 

1 1 

B 2B 2B 2B 


25 

2 


x- (x-9)- (9-x+l)- (27- x+1) 

4 3 323-x 



Bl-x 717 -x 3-x 


3 

3 



#11: chr= 


x .(x+2.^2-4).(x-2-^2-4) 

3 2 

x ■ (1-x)- (x-r£) 

12 -x -(x-6)-^x j-(3-x+l) 


2- (B-B-x)- (4-x -4-X+13) 


2B 


9-x- (x-9) ■ x 


27 


(9-x+l) 


#12: chc= 


x ■ (x+2-^/2-4)- (x-2-^/2-4) 

3 2 

x -(l-x)-(x-rS) 

12 -x -(x-6)-^x j-(3-x+l) 

B ■ (B-B-x) ■ |^x +J3 ■ l j ■ |^x J3 ■ L 

2B 


(9-x+l) 


9-x- (x-9) ■ x+ 


27 


Having some experience in programming it is not really a problem to transfer the idea(s) to 
Tl-Nspire. The screen shot shows the FACTOR-utility. The students have to do factorization 
by hands and then check their results by using the provided functions. 


expand(<ftZ>(3)) 


25*tfrv 2 *z^-400v/7 3 *v 2 *z 2 
240*A 2 vn 4 7? 2 -540*/r 6 vw 2 7? 2 
560-k 7 •n Z -z Z -S7S-hn Z -z 5 


factor 

25*ot*v 2, z^-400*ot 3 *v 2 *z 2 

240*/? 2 v?7 4 ^ 2 -540*/7 6 *ot 2 ^ 2 

560-* 7 V7 3 *z 3 -875*Av7 3 *z 5 


25 , WV 2, Z 2 *(z-2 , /77)*(z+2*/77)*(z 2 +4*W 2 ) 
-60 A 2 ’(s* A 2 -2 m m) m [$ m h ^+2 m m\ m m 2 77 2 
-35*fc/7^*z^*(5*z-4*£ 3 ) , (5 , z+4*/r') 

expand(/?/?()) 

900* ( g 5 *A 3 -900* ( g 4 *A 3 -900* ( g 3 *A 3 

factor(900* < g 5 */f 3 -900-g 4 */7 3 -900* ( g 3 - 

A 3 ) 

900-g 3 -^ 2 -g-l)-A 3 

expand (sA()) 

224-c 8 *^z-28*c 2 */rz 4 

factor) 

i 224*c 8 -^z-28-c 2 -^z 4 ) 

-28*c 2, /r , z*(z-2*c 2 ) , (z 2 +2*c 2 *z+4*c 4 ) 


expand [zu *W) 


"6^g*s 4 -8 ,, i g 3, s 4 -30 , c 2 ’/g*s 3 -40 , c 2 ' i g 3, s 3 

-240*/7^*^ 4 *v 6 -540 , /7 3 *«7 10, v 4 

200* < g 5 -OT 3 ^ 3 +800-,g 3 *OT 3 ^ 5 

8-c Z -f-k-24-c Z -f-k Z +18-c Z -f'k 5 


7/991 


I don't want to demonstrate the Long Division because this can be found in DERIVE News- 
letter 71/72 from 2008 [15] . 
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Many many years ago I wrote a BASIC program for teaching and training set theory. It was 
an “intelligent” program, because I implemented a parser in order to perform any set opera- 
tion with 2, 3 or four sets of numbers or characters (given or randomly chosen). 

|JV|BiChJL J Mi i'iii- i'l QiCi iJ.Ui.ZlJUJ U5.UJ 

3i«GE 94KB Anvmiung 30.10,199002:01 


A N 

G A B E : 



G = 

{9,17,12, 

20,15,2,3,25,7,10,16,11,6, 

,22} 

A = 

{20,16. 

,17 

,11,6,7,15,12,25,22,2,9,3, 

,10} 

B - 

{20,22. 

,6 , 

15,3,9,11,10} 


C = 

{12,25. 

,11 

,16} 



Aufgabe : (A M) \ (A N C) 


(A S B) S (A \ C) = {9,17, 15} 
Eas stimmt leider nicht. 

(A S B) S (A S C) = 116,12,25} 


( E N 1 E R > 


The set operations could be visualized. 


P r q b 1 e 


(AUB) n (AUC) fl (BUC) 



Fur die Mengen nur 
den Buchstaben, fur 
die Teilmengen nur 
die Ziffer eingeben. 

Teilmengen uerden 
griin gefiillt. Beleg- 
te Teilmengen konnen 
durch nochmal iges An- 
sprechen geleert uer- 
den . 

Die Eingabe kann ab- 
geschlossen uerden, 
uenn nur melir dunkle 
Teile vorkommen. 




? lost das Problem ? 


The nice thing is that students in my former school are still using this “antique” program 
which is running in the DOS-environment. 


It was a challenge for me to produce a similar tool with DERIVE - without facing the problem 
to implement again a parser - CAS does the job for me. 


This is the recent version of set theory exercising and learning. The students have the choice 
between basics and more complex operations with 2, 3 or 4 sets. The program provides a 
random generated universal set with the respective subsets and a list of problems. 


Here is the section for tasks involving three subsets: 
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Given is a universal set together with 4 subsets A, B, C and D. 
Find the following sets: - Simplify sols for the solutions! 


(A u B) n (B u C) 

B u (C n D) 

(A u C) n (A u D) n (C u D) , 
(A n B’) n D 
A’ \ (B u C) 

(B \ C) u (B n D) 

C n (A n B)’ 

(A u B)’ n (C n D)’ 


(A n B) u (A n D) 
A n (C u D’) 

(C u D) \ A 
(B \ D) \ (B \ A) 
(A \ B) u (A \ D) 
A' u C' u D' 

(D’ u (A \ C))’ 

(C u B) \ (C \ D) 


#16 : th ree= 


Universal Set: {2 ,S ,11,13 ,16 , 21, 24 , 29 , 30} 

Subset A: {2,3,13,16,21,30} 

Subset B: {8,11,16,21,30} 

Subset C: {2,8,11,13,16,21,24,29,30} 

Subset D: {8,11,16,21,29,30} 


The students shall solve the tasks by paper and pencil and then compare their results with 
the given solutions. Own operations can be entered. It is also possible to generate sets of 
random generated sets for self made exercises. 

I call the solutions for the tasks given above: 


{2,8,11,13,16,21,30} 

{8,11,16,21,29,30} 

{2,8,11,13,16,21,29,30} 


{8,16,21,30} 

{2,8,13,16,21,30} 

{11,24,29} 


#17: sol s= 


{} 

{} 

{8,11,16,21,30} 

{2,11,13,24,29} 

{24} 


{} 

{2,13} 

{2,11,13,24,29} 

{3,11,16,21,29,30} 

{} 


Providing exercises with Venn diagrams was not so easy, but finally I found a way to realize 
this in a satisfying way. 

I show the part with three subsets embedded in a universal set. First see the “instructions” 
together with a list of provided problems. 
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Exercising set diagrams (Venn diagrams) with three sets 
Double click on the provided graph to open the 2D plot window. 

Editing A or B or C and plotting shades the respective sets, editing and plotting ql through q8 
shades the distinct disjoint subsets. 

This might help to visualize the following set operations: 

Sketch the Venn diagrams for the given set operations by paper and pencil - 
- or using the technique described above. 


(1) 

(AuB)nC 

(2) 

(AnC)uB 


(3) 

(AuB)n(Bu C) 

(4) 

(AnB)u(AnC) 


(5) 

An(BuC) 

(6) 

(A u C) n (A u B) n (B u 

C) 

(7) 

(AuC]\B 

(8) 

(A n B) 1 u (A n C)' 


(9) 

A'\(BuC) 

(10) 

(A\B)u (An C) 


(11) 

(A \ C) u (C \ B) 

(12) 

A 1 u B 1 u C 


(13) 

(A u B u C) \ (A n B n C) 

(14) 

(AuB)'n(CnB)' 


(15) 

(A\C)'u(B\C) 

(16) 

((A n B) u (A n C) u (B r 

iC))\(AnBnC) 


The correct answer will be given by pl3(# of the task), eg. pl3(12) plots A' u B' u C. 

You can visualize your own diagrams but you have to take in account that DERIVE cannot plot set operations 
buttheireqivalentas Boolean expressions. It is easy to convert the set operations to Boolean expressions: 


A 

u 

B -► 

A or B 

or A 

v B 

A 

n 

B -► 

A and 

B or A 

A B 

A' 


-e 

not A 

or -.A 


A 

\ 

B -» 

A and 

not B 

or 


Let's assume that student Josef wants to treat problem (10). He marks the disjunct subsets 
q2 and q5 (= A\B) together with q8 (= AnBnC) by entering and plotting [q2,q5,q8], because 
he thinks that this is the solution. Then he plots pl3(10) to check if he is right or not. 



I will skip the Solving Triangles Training (will follow in another DNL) and remind you that the 
utility for working with times was published in DNL #74 (“What's the TIME, Grandie?) 1161 . 

One of my latest products is a tool for analysing conics: 
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#1: LOAD(D:\DOKUS\SCHULE\ski'ns\comcs.rnth) 

#2: start 


Conics Trainer - Trainer fur Kegelschn itte 

coni = Center in origin unrotated - Kegelschnitt in Ursprungs lage 
conZ = shifted origin, unrotated - verschobener Kegelschnitt 
con3 = conic shifted and rotated - Kegelschnitt in allgemeiner Lage 
con = random choice from above - Zufal Isaufgabe aus den obigen 


ans = gives the analysis; ant = liefert die Analyse 


#4: 


#5: 


2 2 
con3 = (20 *x - 38-x-y - 14-x - 14-y - 15-y = 1) 

Type Center Real axis 2a Imaginary axis 2b 


ans 


Hyperbola [-0.06942, -0.4414] 0.7052 


0.6266 


Real Vertices 
0.07450 -0.1195 

-0.2133 -0.7634 


Asymptotes 

1. 372 ■ t - 0.06942 0.6192*t - 0.4414 

0 . 4535 • t - 0.06942 - 1.435«t - 0.4414 


The student performs the calculation - possibly supported by CAS - and checks the ob- 
tained results. As an add on he/she can plot the conic and its center, the vertices and the 
asymptotes in case of a hyperbola. 
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I presented a lot of examples from various fields of school mathematics performed with vari- 
ous pieces of soft- and hardware. So you might have made up your own opinion about Pros 
and Contras for using technology in this very special way. 

What are my personal Pluses and Minuses? I'll start with the MINUSES: 

- Very (too?) close to traditional (fundamental) mathematics education 

- There is a danger that training could become a main activity in math teaching 

I am quite sure that the reader will add some other objections - which I probably will share. 

My PLUSES are: 

+ The tools could (should) be developed, tested and improved in cooperation 
with the students, which asks for deeper understanding of the algorithms and 
techniques 

+ Encourages self responsible exercising parallel to math education (if neces- 
sary) 

+ Hesitating teachers could be led to the use of CAS - and hopefully will not 
use CAS for training manipulation skills only but will later change - or at least 
consider to change - their methods teaching maths. 

+ Students like it! 

Let me add one more aspect which is very important from my point of view: 

The kids/pupils/students have to learn a lot of rules. I believe that nobody will teach calculus 
without talking about the chain rule. The students then know the rule and have in their minds: 
"Outer derivative times inner derivative", which is ok. 

But it is my strong opinion that a special competence is needed to apply this rule in the ap- 
propriate way and to know when just this rule must be applied. And this is valid for many 
cases applying rules - also in daily life. Knowing the rule is not sufficient enough. One must 
be able to recognise when and how to apply the rule - and not in all cases does the com- 
puter all the work. 


Exercised enough?? 
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There is some space left. So I cannot resist to add my Taylor-Series-Exercising Tool: 



The student is given a random generated function together with location xO where to develop 
a Taylor series of requested degree, plot shows the respective graphic representation. 
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Making Life in an Introductory Linear Algebra Course 

Easier with DERIVE 

KARSTEN SCHMIDT 

Schmalkalden University of Applied Sciences, Germany 
email: kschmidt@fh-sm.de 

The Schmalkalden University Faculty of Business and Economics moved its intro- 
ductory linear algebra course from the classroom to the PC lab, and purchased a 
DERIVE license that also allows its use on the students’ own PCs. A DERIVE util- 
ity file was then developed to facilitate exercises with special matrices throughout 
the course. This utility file allows the computation of zero matrices and vectors, 
matrices and vectors of ones, random matrices, as well as idempotent and orthogo- 
nal matrices, “just-in-time” whenever they are needed during the course. This paper 
demonstrates how beneficial it is for students to sit in front of a PC in an introduc- 
tory linear algebra course, from the very basic to the more advanced topics. 


1 Some preliminaries 

A matrix is a rectangular array of elements in m rows and n columns. We denote a matrix by a bold- 
face capital letter, e.g. A. In DERIVE , a matrix can be defined simply by clicking on H, setting the 
number of rows and columns, and entering the elements. An element of a matrix is denoted by the 
corresponding lowercase letter with a double index, the first index being the row, and the second the 
column index. For example, the element in the second row and first column of A is denoted by a 21 . 

A matrix with only one column ( n = 1 ) is a (column) vector. Hence, throughout the course, a vector 
has to be defined in DERIVE as a matrix with one column; the vector data type (e.g. defined by click- 
ing on @) is not used. We denote a vector by a boldface lowercase letter, e.g. x. Row vectors are not 
defined; transposed column vectors are used instead. A matrix is transposed by taking its rows (col- 
umns) and writing these as columns (rows) of the transposed matrix. The transpose of a matrix or vec- 
tor is denoted by a prime, e.g. A ' . Clearly, the transpose of an mxn matrix is an nxm matrix. In 

DERIVE the transposition sign is a grave accent which can be entered, for example, by clicking on I 

in the symbol list at the bottom right of the DERIVE window. 

A matrix which has the same number of rows and columns ( m = n ) is called a square matrix. The 
diagonal from the upper left to the lower right of a square matrix is called the main diagonal. 
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2 Matrices and Vectors Containing Os and Is 

It is important for students to develop a good understanding of 

• matrices and vectors containing only Os (zero matrices and vectors), 

• matrices and vectors containing only Is (matrices and vectors of ones), and 

• square matrices containing Is on the main diagonal and otherwise Os (identity matrices). 

While zero matrices and vectors behave similarly to the number 0 e □ , matrices and vectors of ones 
do not behave similarly to the number led, but rather identity matrices are the matrix algebra ana- 
logues to 1 e □ . However, it is rather tedious to verify this in exercises. The utility file therefore in- 


cludes the following five functions: 

0(m , n) generates an m x n zero matrix (1) 

o(n) generates an n x 1 zero vector (2) 

J (m , n) generates an m x n matrix of ones (3) 

1 (n) generates an nxl vector of ones (4) 

I(n) generates an nxn identity matrix (5) 


The names of the functions were chosen such that they are close to commonly used symbols for such 
matrices: O for zero matrices, J for matrices of ones, I for identity matrices, o for zero vectors, and “1” 
(lowercase L) for vectors of ones (since this resembles the symbol 1 which is typically used for vectors 
of ones). It is important for the students to understand that if they use the utility file (i.e. if the file is 
loaded into DERIVE ), they cannot use the above names for anything else (e.g. for a matrix), as this 
would overwrite the function in the utility file. Note that the I (n) -function was added simply to have 
a shortcut notation for DERIVE’ s built-in function for the computation of identity matrices, 
IDENTITY_MATRIX(n). 

The utility file also contains a function that generates an m x n random matrix of nonnegative integers 
(cf. next section). The third parameter of this function determines the highest possible integer such that 
the following function is actually within the scope of this section: 

RNDU (m , n , 1) generates an m x n random matrix of zeros and ones (6a) 

We define the matrix 

^12 3 

A = 4 5 6 

3x3 

^7 8 9 

for which we want to check some of the following properties of zero matrices: 

A 0 = 0 ; O A = O ; O = o o' ; 0 0 = 0. 

mxnnxl mxl Ixmmxn Ixn mxn mxllxn mxnnxp mxp 
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Note that the last of the above properties implies that the product of any square zero matrix with itself 
is a square zero matrix of the same dimension, making square zero matrices a standard example for 
idempotent matrices (any square matrix A with AA = A is called an idempotent matrix, cf. section 4). 
Screenshot 2.1 shows a few exercises with zero matrices. 

#1: LOAD £C:\Prog ram Files\TI Education\Derive 6\Math\SpeM.mth) 

12 3 


#2: A := 


4 5 6 

7 S 9 


#3: 


#4: 


#5 : 


#6: 


A -0(3 r 2) = 


0 0 
0 0 
0 0 

□(3)’ -A = [[0 P 0 r 0]] 
0 0 0 


□(2) = 


0 ( 2 - 2 ) - 0(2 r 2 ) = 


0 0 0 
0 0 
0 0 


Screenshot 2.1 


We continue with matrices and vectors of ones. For the previously defined matrix A we want to check 
some of the following properties: 


1 xn nx 1 


i = 1 


n \ 

IX 

7=1 


a' 1 = 1 'a = ; 


A 1 = 

mxn nx 1 


z 

V7'=! J 


a mj 


r a = 

1 xm mxn 


IX ••• 

V i= 1 i=l J 


'l: 1 ' 


1 xm mxn nx 1 


»=1 7=1 


7 = 11 '; 

mxn mx 1 1 xn 


J J =n J 

mxnnxp mx p 


Clearly, the first four properties show that vectors of ones are useful for the computation of sums, 
namely of the elements of a vector, or the row, or column, sums of a matrix, or the sum of all elements 
of a matrix. Screenshot 2.2 shows a few exercises with matrices of ones. 



Screenshot 2.2 


For a change, the summation of the elements of a vector a is done for a 2x1 vector of arbitrary ele- 
ments, which are denoted by al and a2. The result of the multiplication of two matrices of ones is 
unexpected by many students, as it is not another matrix of ones. 

After realizing that matrices of ones do not behave similarly to the number led, attention can im- 
mediately be directed to identity matrices. For the matrix A and the vector a already defined we want 
to check the following property (see Screenshot 2.3): 


A I = I A =A . 

mxn nxn mxm mxn 



It is also shown that the product of the 2x2 identity matrix with itself is the 2x2 identity matrix. As 
this is true for any identity matrix (which follows immediately from the above-mentioned property by 
choosing A = I ), identity matrices are another standard example for idempotent matrices. 
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3 Random Matrices 


Many students appreciate having the option to use additional example matrices to practise certain top- 
ics in matrix algebra. The utility file therefore includes the function 

RNDU (m , n , max) generates an mxn matrix of random nonnegative integers (6) 

The random numbers are from the set {0, 1, 2, . . . , max } , i.e. the third parameter determines the highest 


possible value; for example, if max = 1 , the generated matrix contains only Os and Is (cf. Screenshot 
3.1). The RNDU-function should only be called with a positive integer as the third parameter. If 
max < 0 , an error message is printed, if max > 0 but not an integer, no random matrix is generated. 
Note that the RNDU-function uses the built-in DERIVE function RANDOM(n) with n = max + 1 . 


# 2 : 


#3: 


#4: 


RNDU(2 ? 2 r 1) 
RNDU(2 ? 2 , 1) 
RNDU(2 F 2 p 1) 


0 1 ' 
0 1 . 
1 1 ' 
1 lj 
1 1 ' 
1 0 J 


Screenshot 3.1 


Only the last of the three matrices generated in Screenshot 3.1 has full rank (is nonsingular), the other 
two have rank 1. As any element of a matrix generated by the RNDU (2 , 2 , 1) -function is either 0 or 1, 
and as there are only 16 different 2x2 matrices with two different elements, the probability of getting 
a matrix having full rank is lower (% 6 ) than it is of getting a singular matrix ( % ). This allows an 
interesting discussion of the rank of a matrix. 

In general, the probability of getting a rank-deficient matrix decreases if the number of rows, the num- 
ber of columns, or the max value is increased. For example, all three matrices in Screenshot 3.2 are of 
full row rank. 


#5 : 


#6: 


#7: 


RNDU(2 f 3 , 9) 
RNDU(2 ? 3 P 9) 
RNDU(2 f 3 r 9) 


7 5 9 

0 9 6. 

2 0 6 ' 
17 8 . 

4 2 1' 

9 3 0. 


Screenshot 3.2 


Note that there is also a built-in DERIVE function similar to the RNDU-function: the 
RANDOM_MATRIX(m, n , s) function generates an mxn matrix of random integers from the set 

s + 1, —s + 2, . . . , — 1, 0, 1, . . . , s — 1} . 


p 44 


Karsten Schmidt: Introductory Linear Algebra Course 


D— N— L#76 


4 Idempotent Matrices 


As we have now already come across two examples for idempotent matrices (any square matrix A with 
the property AA = A), namely square zero matrices and identity matrices, we should spend a little 
more time with these matrices. The utility file includes the function 

CNTRC n) generates an idempotent nxn matrix (7) 

The formula used in this function, 

nxn nxn 


generates a so-called centering matrix (using an identity matrix and a square matrix of ones). Center- 
ing matrices are used, for example, in multivariate statistics. Screenshots 4.1 and 4.2 show what cen- 
tering matrices look like and also that the 2x2 and 3x3 centering matrices are indeed idempotent. 


# 2 : 


#3: 


1 1 


CNTR(2) = 


2 2 
1 1 


2 2 


CNTR(2) - CNTR(2) -CNTR(2) = 


0 0 
0 0 


Screenshot 4.1 



Obviously, the CNTR-function can only generate one particular idempotent matrix of a specific dimen- 
sion. Therefore, the utility file contains another function for the computation of an idempotent matrix, 
the IDEM-function, which uses a property that holds for the Moore-Penrose inverse A + of a matrix 

nxm 

A : 

mxn 


A + A and AA + (as well as I-A + A and I-AA + ) are idempotent matrices. 
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The Moore-Penrose inverse of any matrix A is the (unique) matrix A + satisfying the four condi- 

mxn nxm 

tions 


AA + A = A ; A + AA + =A + ; (a + a)' = A + A; (aA + )' =AA + . 


As DERIVE does not include a function for the computation of the Moore-Penrose inverse, the utility 
file includes the two functions 


MPIV(a) 

computes the Moore-Penrose inverse of any n x 1 vector a 

(8) 


MPI(A) 

computes the Moore-Penrose inverse of any mxn matrix A 

(9) 


Both functions are described in detail in Schmidt (2003). 

If A is a square non-singular matrix, the Moore-Penrose inverse and the inverse A -1 coincide, and we 
have A + A = A -1 A = I and I - A + A = I - A~ l A = O , i.e. we are back to the two standard examples 
for idempotent matrices from section 2. Therefore, the more interesting cases are when A is either 
square but singular, or non-square. The matrix A from section 2 is such a singular matrix. The Moore- 
Penrose inverse of A is computed in Screenshot 4.3. 




r i 

2 

3 ' 





#6: 

A := 

4 

5 

6 







. 7 

S 

9 . 


23 

1 

11 ' 








36 

6 

36 


#7: 





MPI(A) = 

1 

IS 

0 

1 

IS 








19 

1 

7 








36 

6 

36 . 



Screenshot 4.3 


The utility file includes the function 


IDEM (A) generates an idempotent nxn matrix from any mxn matrix A (10) 


5 

1 

1 

6 

3 

6 

1 

1 

1 

3 

3 

3 

l — 1 

1 

5 

i 

cr> | 

3 

6 


#S: 


I DEM (A) = 


# 9 : 


IDEM(] (3 f 2)) 


Screenshot 4.4 
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The formula used in this function is A + A . Screenshot 4.4 shows the idempotent matrix which is gen- 
erated when matrix A is passed as the parameter, and also the one generated by passing the 3x2 ma- 
trix of ones. 


Note that finding the rank of an idempotent matrix is a relatively easy task, as rank and trace are iden- 
tical in this case; for the first example in Screenshot 4.4 we get | + f = 2 , for the second \ + \ — 1 . 
We can also quickly double-check this on the PC (Screenshot 4.5): 


#10: 

RANK (IDEM (A)) = 2 

#11: 

TRACE (IDEM(A)) = 2 

#12: 

RANK(IDEMCJ(3, 2») = 1 

#13: 

TRACE (IDEM(I (3 , 2)» =1 


Screenshot 4.5 


5 Orthogonal Matrices 


Finally, we want to consider orthogonal matrices (any square matrix A with the property A! — A * 1 ). 
The utility file includes the function 


ORTH (a) generates an orthogonal n x n matrix from any n x 1 vector a ^ o (11) 


The formula used in this function, 

I -2 a a' , 

nxn nx 1 1 xn 


# 2 : 


b :z 


5 

4 


#3: 


#4: 


b’-b = [[1]] 


7 24 


ORTH(b) = 


25 

24 


25 

7 


25 25 


Screenshot 5.1 


generates an orthogonal matrix if a' a = 1 . Hence, the vector a which is passed as the parameter will be 
transformed within the ORTH-function such that it is of length 1. Let us nevertheless start with a vector 
b that has length 1 anyway, as in Screenshot 5.1. 
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In order to check if this matrix is indeed orthogonal (cf. Screenshot 5.2), we have to compute both its 
inverse and transpose, and see if these two matrices are identical. Since for any orthogonal matrix we 
have 


A' A = A~ l A = I , 


a second method to prove that a matrix A is orthogonal is to show that A' A = I . 


#5 : 


#6: 


-i r o o 

ORTH(b) - QRTH(b) 1 = 

L o o 


ORTH(b) 1 -QRTH(b) 


" 1 0 

. 0 1 


Screenshot 5.2 

For a second example (Screenshot 5.3) we choose the 2x1 vector of ones as the parameter of the 
ORTH-function. 


#7: 


#S: 


#9: 


ORTHO (2 » = 
-1 


0 -1 

-1 0 


QRTH(1(2)) - QRTH(1(2)) 1 

ORTH (1(2)) 1 -ORTH (1(2)) = 



0 0 

0 0 

0 ' 
1 . 


Screenshot 5.3 


To complete this section we take advantage of the graphical capabilities of DERIVE. Any column of 
an orthogonal matrix has length 1, and any two are pairwise orthogonal. Since both examples are two- 
dimensional, its column vectors lie on the unit circle and form a right angle (Screenshot 5.4). The 
graph also shows the vector b from the first example. 



Screenshot 5.4 
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Remarks 

This paper is linked to a presentation given on June 26, 2009, at the ACA 2009 Conference in Mont- 
real, Canada. 

Some portions of this paper are also published in the Proceedings of the 7 th Delta Conference in 
Gordon’s Bay, South Africa. 
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